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Evolute 

Involute 

Wavefront 

Caustic surface of two branches (co,, co 2 ) 

Incident and reflected rays 
Normal to a surface 
Lines of curvature 

Characteristic directions (eigenvectors on a wavefront) 

Characteristic values (eigenvalues), principal normal curvatures 

Normal, Gaussian, and mean curvatures 

Generic parametric surface coordinates 

Tangents to a surface 

Tangent rates 

Coefficients of the First Fundamental Form 
Coefficients of the Second Fundamental Form 
Coefficients associated with k„ e, 

Azimuthal angle in Euler's Theorem 

Derivatives in Rodriques’ formula 

Principal normal radii of curvature (Flux-density context) 

Flux-density 

Angles of incidence on deflector and receiver 

Coefficients of the Flux-density formula 

Distance from reflector to observer (Flux-density formula) 

Focal length of a paraboloid 
Reflector surface 

Wavefront aberration along normal to Gaussian sphere (GS) 
Reflector distortion 
Optical path length 

Planes associated with a Darboux frame 

Centers of curvature 

Edges of regression (space curves) 

Aberration coefficients 

Optical invariants (h 2 , a 2 , h o cos f) 

Azimuthal variable 
Radial variable 
Radius of Gaussian sphere 
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ANALYTICAL CAUSTIC SURFACES 


INTRODUCTION 

The study of caustic surfaces is rooted in antiquity. “One of the earliest discoveries in optics (F. Maurolycus, 1575) 
was that the rays of a normal system are tangential to a surface, the so-called caustic surface, and that the general ray 
is tangential to the caustic surface in two points,” Ref. 1, p. 156. The study of caustic surfaces persists to the present. 
"Thom himself observed that a usual interpretation of his theorem on the classification of singularities of gradient 
mappings as elementary catastrophes were the so-called caustics.” Ref. 2 (1982), p. 1080. A great amount of material 
on caustics can be found between the discoveries of Maurolycus and the observation of Thom. Many statements per- 
taining to caustics stir the imagination; some are provocative. For example, the caustic has been described as “a zero, 
one-, or two-dimensional manifold in every neighborhood of which more than one ray of any infinitesimal tube of the 
field passes through some point,” Ref. 3, p. 326. In geometrical optics: “The envelope to a family of rays reflected 
from some surface is called a caustic, and the Hilbert integral along an envelope degenerates into a variational in- 
tegral,” Ref. 4, p. 79. From differential geometry: “The center surfaces are the loci of the centers of principal (normal) 
curvature on the surface,” Ref. 5, p. 94. The caustic surfaces are usually associated with focusing, Ref. 6, p. 259. 

The objectives of this document are limited to a discussion of several alternative means for obtaining analytical 
caustic surfaces and the exploration of wavefronts and caustics of a general form that derives in a natural way via 
aberration theory. An overall appreciation of the structure of this document may be obtained from Figure 1 . A family 
of lines defining the caustic surfaces is obtained analytically from rays via the law of reflection, perhaps the simplest 
and most familiar approach. Introduction of the optical path length, in conjunction with the law of reflection, leads to 
the wavefront approach. Procedures based on differential geometry permit the analytical development of center sur- 
faces, or caustic surfaces, as a set of points via wavefront normals and principal normal curvatures (k,, k a ). A flux- 
density algorithm obtains an equation for the caustic surfaces in terms of the intrinsic Gaussian, mean, and normal 
curvatures (Kg, Km, Kn) of a reflector by exploiting a singularity of the formulation. This approach allows the mapp- 
ing of contours of equal flux density and the imbedding of the caustic surfaces. A second method that leads to 
parametric equations of the caustic surfaces uses a tensor algorithm that regards only the source and the reflector, to 
the exclusion of the wavefront. 

The practical or applications aspect of “Analytical Caustic Surfaces” is implicitly contained in Figure 1 . It can be 
seen that provision is made for a distortion of the reflector, resulting in a distorted wavefront. Alternatively, an 
idealized wavefront may be assumed at the outset, and caustic surfaces may be obtained after superimposing combina- 
tions of many types of appropriately weighted classical wavefront aberrations. Sources (S,) may be generalized to 
represent converging, diverging, or planar types, and superposition of sources is allowed under the flux density 
algorithm. Intermediate media, between source and reflector, may be treated as a modification of the illumination 
block of Figure 1. 

The relationship between rays, wavefronts, and caustics is illustrated for convenience in two dimensions in Figure 
2 and initial discussion is in terms of involutes and evolutes. Actually, the latter are defined for space curves in three 
dimensions as follows: 

(1) The evolute of a curve is the geometric locus of the centers of curvature of a curve. Ref. 7, p. 58 
and 

(2) The involute of a curve lies in the tangent surface of that curve and intersects the tangent lines orthogonally. Ref. 

8, p. 84. 

Rays are seen to be orthogonal to wavefronts and tangent to the caustic curve. The arc distance S« on the caustic 
equals the tangent distance from the wavefront to the caustic. It is therefore the inverse of the curvature at every point 
of the wavefront. The evolute of the wavefront is the caustic. The involute of the caustic is the wavefront. The 
wavefront is not unique and depends on the choice of an initial point of the caustic. Bertrand curves are recalled. Ref. 
8, p. 98 


EV(A) = fi 

(1) 

INV(O) = A 

(2) 


Generation of Figure 2 is made easy since the caustic (Q) was taken to be a circle for which a closed-form expres- 
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Figure 1. Analytical Caustic Surfaces. 
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Figure 2. Rays, wavefronts and caustics. 



sion of arc length is available. This in turn provides the radii of curvature of A directly. The example is a particularly 
good one since a portion of the Tschirnhausen cubic is nearly circular. The latter curve is a cross-section of the caustic 
surface for off-axis plane-wave incidence on a parabolic arc. Background material on the derivation of this cubic 
curve is found in Ref. 9 and Ref. 10. Figure 3 illustrates the formation of Tschirnhausen’s cubic by tangent rays and is 
a point of departure for developing analytical three-dimensional caustic surfaces of paraboloids and other reflectors 
under various assumptions regarding illumination and distortion.' 

The objectives of this document are brought out in greater detail in Figure 4. Representative wavefronts (A) and 
their associated center or caustic surfaces (w„ « 2 ) are shown without any reflector surface (y) or illuminating rays (S,). 
Since the system of the reflected rays (S r ) is orthogonal to the wavefront (A), the wavefront normals and the reflected 
rays are co-directional or at most anti-directional. The rays (S r ) are tangent to the caustic branches (w,, o> 2 ). An 
analytical means of determining the caustic surfaces may therefore be reduced to the problem of finding the wavefront 
normals (N), a family of lines in three-dimensional space. Alternatively, the caustics may be found as a set of points, 
using the definition of a center surface, if the extrema (k„ k 2 ) of normal curvatures (k) can be found at every point P 
on the wavefront A. Standard methods of differential geometry are employed to obtain both the unit normals (N) to a 
surface (x), via the un-normalized tangents (x„, x„), and the principal normal curvatures (k,, k 2 ). These methods also 
permit the writing of equations, practical closed-form expressions, to describe the caustic surfaces analytically in 
many cases. An advantage of obtaining the caustic surface equations is that they may in turn be analyzed for cur- 
vature. 

An objective of this document is to include wavefront aberrations and pursue the analysis when the tangents x u and 
x, are no longer orthogonal, and all terms of the first and second fundamental forms of differential geometry are dif- 
ferent from zero. The directions of the principal lines of curvature (LC„ LC 2 ) imply orthogonal vectors (eigenvectors) 
associated with the principal normal curvatures (eigenvalues) and principal centers of curvature (Ci, C 2 ) in Figure 4. 
The vectors N, e ( , e 2 then constitute an orthogonal moving triad on the wavefront (A), and the pairs of vectors (N, e 2 ), 
(N, e 2 ) identify with (principal) planes of curvature. 1 2 All of the preceding discussion pertains to the formation of 
caustic surfaces from a wavefront using standard methods of differential geometry to obtain the characteristic values 
(ki, k 2 ) and the normal (N) of the latter surface. 

Similar techniques may be applied directly to a reflector surface instead of a wavefront to obtain the caustic sur- 
faces. Characteristic values and surface normals are required as before. The stipulation of an incident ray (S,) and the 
reflector normal (N) defines a plane of curvature intersecting the reflector and therefore a normal curvature (Kn), not 
necessarily a principal normal curvature. The characteristic values (k„ k 2 ) of the reflector determine the Gaussian (Kg) 
and mean (Km) curvatures. An existing analysis, predicated on these intrinsic curvatures (Kg, Km, Kn) of the reflector, 
provides both flux density and caustic information through exploitation of a singularity of the formulation, as men- 
tioned, previously. Standard methods of differential geometry may therefore be used in conjunction with alternative 
theories for obtaining caustic surfaces. 


1 Once the nodal cubic has been identified it may be rewritten in parametric form as 

x = 3a(t*-3) , y * at(t*~3) , 
from which the expression for arc-length, 

S = at(t 2 +3) 
follows. 

2 The triad is called a Darboux frame for surfaces, and is the analogue of Fresnel’s frame (t, n, b) for space curves. Ref. 14, p. 210; Ref. 5, p. 261 . 
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Figure 3. Tschirnhausen’s Cubic (caustic arc) for 0- = 10°. 
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Figure 4. Center Surfaces (ojj, oo 2 ) of a wavefront (A). 


6 



CENTER SURFACES 

A review of the theory of center surfaces is useful before setting out to find the surface normals and principal nor- 
mal curvatures of either wavefronts or reflectors. From the three fundamental equations of differential geometry, 


I = dx • dx (3) 

II = -dx • dN (4) 

III = dN • dN , (5) 

the normal curavature is defined as 

_ II _ — dx • dN _ edu 2 + 2f du dv + gdv 2 (6) 


I dx • dx Edu 2 + 2F du dv 4- Gdv 2 , 


where the surface is: 


and the surface normal is 


x = x(u, v) 


N = N(u, v) 


(7) 


( 8 ) 


Results obtained for principal normal curvature may be verified using Rodriques’ formula since it characterizes the 
eigenvectors. 


dN = -k dx 


In the preceding. 


and 


Here 


and 


since 


dx = Xu du + Xv dv 

dN = Nu du + Nv dv 

E = Xu • Xu, F = Xu • Xv, G = Xv • Xv 

e = L = -(Xu • Nu) = (Xuu Xu Xv)/(EG - F 2 ) ,/2 
f = M = - (Xu • Nv + Xv . Nu) = (Xuv Xu Xv)/(EG - F 2 ) ,/2 
g = N = -(Xv • Nv) = (Xvv Xu Xv)/(EG - F 2 ) 1/2 

Xu • N = 0 , Xv • N = 0 


(9) 

( 10 ) 

(ID 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
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and 


xu • Nv = Xv • Nu 


(17) 


The (E, F, G) and (e, f, g) are coefficients of the first and second fundamental forms. 

The unit surface normals are obtained from a cross-product of tangents, not necessarily orthogonal, as indicated 
in Figure 4. 


N = (Xu xXv)/ 1 Xu xXv | = (Xu xXv)/(EG • Fa) 1/2 


(18) 


Letting 


X = dv/du 

the normal curvature of a surface becomes 


e + 2 fX 4- g\ 2 
E + 2FX + GX 2 


(19) 


( 20 ) 


The extrema of k are characterized by 

dk/dX = 0 


( 21 ) 


so that 

II _ f + gX, _ e + fX, . . _ , „ (22) 

K* — — >1 Ij * 

I F + GX, E + FX, 

See Ref. 5, p. 80; Ref. 21, p. 287. 

The X, may be found from 

X, + X 2 = -(Eg - Ge)/(Fg - Gf) (23) 

X, X 2 = (Ef - Fe)/(Fg - Gf) (24) 

See Ref. 4, p. 142 

The principal directions (eigenvectors) associated with the principal normal curvatures (eigenvalues or 
characteristic values) are given by 

e, = (Xu + Xv X,)/H, , i = 1 ,2 (25) 


where 


H, 2 = E + 2F X, + GX, 2 
and are always orthogonal since 

e, • e 2 = (Xu + Xv X.) • (Xu + Xv Xa)/H,H 2 

= [E + (X, + Xa)F + XjXaG]/HiHa = 0 


(26) 


(27) 
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The normal curvature in an arbitrary direction takes the form 


k = k, cos 2 a + k 2 sin 2 a 


and is known as Euler’s tfieorem. 

The derivatives Nu, Nv appearing in Rodriques’ formula are given by 


Nu = 


fF - eG 
EG - F 2 


Xu + 


eF - fE 
EG - F 2 


Xv 


and 


Nv = 


gF - fG 
EG - F 2 


Xu + 


fF-gE 
EG - F 2 


Xv 


( 28 ) 


(29) 


(30) 


and are known as the Weingarten equations. Ref. 5, p. 108. 

A normal curvature (k) is a principal normal curvature (ki, k 2 ) if and only if k is a solution of the quadratic 

(EG - F 2 ) k 2 - (Eg + Ge - 2Ff) k + (eg - f 2 ) = 0 (31) 


Ref. 8, p. 183. 

Since the roots of this equation are 

k (Eg + Ge - 2Ff) * V(Eg + Ge - 2Ff) 2 - 4(EG - F 2 )(eg^T) , (32) 

' 2(EG - F 2 ) 


the mean curvature is 


Eg + Ge — 2Ff _ ki + k 2 
Km = 1 „ " — - 

2(EG - F 2 ) 2 

and the Gaussian curvature is 

Ko = 68 ~ ^ = k, k 2 
EG - F 2 

so that Equation (28) may be rewritten as 

k 2 - 2 Km k + Kc = 0 


(33) 


(34) 


(35) 


The surface normals along a line of curvature are tangent to two space curves, or edges of regression, at distances 
k, _1 and kj*’ from a point (P) of a given surface. Ref. 5, p. 94; Ref. 21, p. 25. Lines of curvature are always or- 
thogonal. Center surfaces (caustics of two sheets) are the loci of the centers of principal normal curvature of a given 
surface (wavefront). Figure 4. 

The points on a surface are classified as follows: 

Elliptic: eg - f 2 > 0 — Kg > 0 (36) 

Hyperbolic: eg - f 2 < 0 — Kg < 0 (37) 
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Parabolic: 

Planar: 

Umbilical: 


Ref. 8, p. 177- 


eg — f 2 = 0 and e 2 + f 2 + g 2 ^ 0 — Kg = 0 

e = g = f = 0 — Kg = 0 

e f g 

k = — = — = — = constant for all directions 
E F G 

p. 184; Ref. 16, p. 96; Ref. 6, p. 261 


(38) 

(39) 

(40) 
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CAUSTICS FROM DISTORTED WAVEFRONTS 


The computational route for obtaining caustics as a set of points from a distorted wavefront is now presented in 
some detail. A cylindrical net or (a, f) parametrization is used to describe an idealized spherical wavefront. Aberra- 
tions are then introduced along the normal of a Gaussian wavefront and, alternatively, parallel to the system axis 
when appropriate to achieve simplification of the expressions. Since caustic information is to be obtained as a set of 
points belonging to center surfaces, machine computations are anticipated at an early stage. The partial derivatives 
Xu, Xv, Xuu, Xuv, Xvv are lengthy, but usually manageable algebraic forms. Forming E, F, G analytically involves 
vector dot-products. Forming e, f, g involves both vector dot- and cross-products. The latter are not developed 
analytically. 

A physical picture associating reflector, wavefronts, and caustic (u t ) is given as Figure 5. The wavefronts (A,) cor- 
respond to reception of a 10-degree plane wave by a paraboloidal reflector, and are illustrated in two dimensions only. 
Caustic analysis will, however, be in three dimensions throughout this document. Figures 2 and 3 are closely related to 
Figure 5. The cusps of the wavefronts appear to lie on the caustic; this observation is not explored here. Direct deter- 
mination of the center surfaces or caustics of wavefronts such as Ai, A», A j is sought for a large class of problems. 
Figure 5 represents a specific choice of source and surface and is included here for illustrative purposes only. 

Introduction of the Gauss/Seidel aberrations of degree 2/4, to the exclusion of the Schwarzschild aberrations of 
degree 6 and other aberrations, is arbitrary. The method of center surface determination is not restricted to these 
classic aberrations; they simply provide a useful set of examples encountered in optics and microwaves. The Oaussian 
sphere may be associated with the class of problems that regards circular apertures. It is replaced, for other 
geometries, by conformable reference wavefronts. Cylindrical and toroidal geometries are examples where the Gaus- 
sian sphere is inappropriate. Ref. 17, Ref. 18. 

Appendix A and Appendix B contain the details pertaining to distorted wavefronts for the Gauss/Seidel aberra- 
tions when perturbations are directed parallel to the system axis (Z) or along the normal (N*) to the Gaussian sphere, 
respectively. The latter case is illustrated in Figure 6. It can be seen that the set of normals for f equals a constant on 
the Gaussian sphere is coplanar. This is not necessarily true for the corresponding set of normals on the distorted 
wavefront. As a result, a planar arc on the Gaussian sphere maps to a pair of three dimensional space curves or edges 
of regression lying on the caustic surface branches (to,, w 2 ). A caustic surface may be represented by a sufficiently large 
number of such curves (T) in (x, y, z) space. A single curve may be dealt with by considering projections onto the plane 
f equals a constant and a plane orthogonal to it by forming 

Ri,2 iNa • k = Ri , 2 ,N a • (0, 0, 1) , (41) 

R „2 ,Na«,N, = R „2 ,Na •(- sin f, cos f, 0) , (42) 

R.,i ,Na • ,Nl = R .,2 ,N a • (cos f, sin f, 0) , (43) 

where ,N a , ,N | , ,Nl, are unit normals to the distorted wavefront (A) and the projection planes, respectively. Here 
R, , 2 are the principal normal radii of curvature. 

Three dimensional isometric representation of the caustic surface branches (w,, o> 2 ) is to be preferred to the projec- 
tion described above. 
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Figure 5. Reflector ( 7 ), Wavefronts (Aj), and Caustic (ojj). 
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Figure 6. Mapping of Surfaces GS: A: ft (cjj, 0 l> 2 ). 
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CAUSTICS FROM FLUX DENSITY SINGULARITIES 


Caustic surfaces may be Obtained from an algorithm that permits the computation Of flux density at generic field 
points with the exception of the caustic coordinates. Refs. 11, 12, 13. Informatidh regarding the caustics is obtained 
by exploiting a singularity of the flux-density formulation 

„ _ C„,r 1 COS <p, COS 4>, . (44) 

rtsi-tn — ' 

|a 0 4 flirj 1 + a 2 r a 2 1 

The denominator of equation (44) is set equal to zero, and the resulting quadratic is solved for distances lta, Ra along 
a reflected ray. Wavefronts do not enter explicitly, and caustics are obtained under this approach in terms of the in- 
trinsic Gaussian, mean, and normal curvatures (Ko, Km, Kn) of a deflector surface and a specified source of illumina- 
tion. Kn is determined by the intersection of the plane formed by the deflector normal and the incident ray with the 
deflector surface. 

In the preceding equation C«, r is a reflection or transmission coefficient associated with a deflector, I is the irra- 
diance due to a source, <p, and \j>, are angles of incidence on the deflector and receiver, respectively. For the present 
discussion, reflectors only, 


C»,r * C* * 1 .0 . ( 45 ) 

The coefficients a 0 , a,, a 2 for reflection are 

d 0 = cds <Pi (46) 

a, » -2(2 Km cos 1 <p, 4 Kn sin 2 <fii) 4- 2 cos v>/(r t ) (47) 

flj =* 4 Kg cos <fi, -2(2Km cos 1 <pi + Kn sin 1 <p)/(x x ) + cos ^./(r,) 2 (48) 


Several specializations are made for the discussion that follows. Incidence on a planar receiver (film) is not an- 
ticipated, therefore the factor cos is set to unity to obtain the flux-density field. Reflection is assumed to be perfect 
and irradiance is arbitrarily taken to be unity also. For illumination by a plane wave, the distance (r,) from source to 
reflector is unbounded so that the following simplification is achieved: 


« _ cos 

FdSi-dSl — 11 “* — “ 1 * “ 

| cos <p x -2(2 Km cos 1 <p, 4 Kn sin 1 y> ( ) r 4 4 4(Ko cos y>,)r 2 2 j 


Ti = R 21 , R 22 


(2 Km cos 2 <p, 4 Kn sin 2 ± V(2Km cos 2 V, 4 Kn sin 2 <PtY - 4Ko cos 1 <p, 

4Kg cos ip. 


(49) 

(50) 


The utility of the flux-density method may be demonstrated by several examples. Consider the paraboloid and, 
following the sign conventions of Ref. 1 1 , take the principal normal curvatures of the refelctor surface to be 


k, 


— 4F 1 

(4F 2 4 a 2 ) 2 ' 2 


(51) 


and 


(4F 2 4 a 2 ) 1 ' 1 

using standard procedures of differential geometry. Ref. 17, Ref. 5. 


( 52 ) 
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When a plane-wave is axially directed upon a paraboloid rotational symmetry exists and 


Kn = k, 


( 53 ) 


for every planar cut containing the system axis. From this 

IC... ~~ 1 (8F + ^ ’ (S4) 

2 (4F 2 + o’) 3 ' 2 2 


and 


Kp = 


4F* 

(4F 2 + e 2 ) 2 


= ki kj 


(55) 


It can be shown that 


, 4F 2 

COS J <Pi = m ' ■ 

4F 2 + a 2 

(56) 

and 


j (T 2 

4F 2 + a 2 

(57) 

so that the radical of equation (50) vanishes. The roots of r 2 are then equal. 


. . , , (4F 2 + a 2 ) 

|r,| - |R»,| - |R»| - 4p 

(58) 

Recalling the equations of a parabola 


z =— -F 
4F 

(59) 

P ** 2F /(I + cos 2 v>() i 

(60) 


and 


cos 2 <Pi = 


N 

(a 2 + z 2 ) ,/2 


(61) 


the result 


P = |r 2 | (62) 

can be shown. The interpretation is that all |r»|, directed along a reflected ray, terminate on the focal point of the 
paraboloid. That is, the caustic surfaces w, and w 2 degenerate to point caustics that are congruent with the focus (F), in 
the present analysis, for an axially-directed plane-wave source. See Ref. 6, p. 259. 
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When the incident plane wave is not axially directed, equation (50) is evaluated. For every value (a) of the radial 
net the local incidence angle (<p t ) at 7 ma Y be found using 


(1N7) • ( - A) = ,N T • S, = cos (63) 

(concave side) (convex side) 

where both 1N7 and iS, are unit vectors. 

An xz-plane cut through the caustic surfaces co 2 for a 10-degree plane wave source is shown in Figure 7. Tschir- 
nhau sen’s cubic (a>i)_is recovered directly, without any representation of the received wavefront as a set of points along 
the reflected rays (S r ). A second caustic branch (o> 2 ) is also obtained here, and Figure 7 should be compared with 
Figure 3. The set of reflected rays in the xz plane of Figure 3 did not provide any information about caustic branch 
(oo 2 ); the three dimensional representation is required. __ _ 

Figure 7, for the paraboloid, may be obtained by plotting R21 x S r and R22 iS r from reflector (7) to caustic sheets (coi, 
u) 2 ). The R2J, R22 are effectively eigenvalues of latent wavefronts which intersect (7) at every (a, f) pair. See Figure 5. A 
more expedient approach is to plot vectors 

x(«.) = x(P) - (R21) iS r (64) 

and 

x(a> 2 ) =x(P) (R22) ]S r (65) 

from the origin (0). Similarly, the coordinates of a computed flux-density (FD) may be written as 

x (FD) = x(P) - (R 2 ) jS r . (66) 

If R 2 is adequately subdivided (AR 2 small), the contour mapping problem of the authors Burkhard and Shealy may be 
relieved by forming a “look-up table” for each ray ( 2 S r ) so that preselected values may be located directly. This ob- 
viates the need for interpolation estimates minimizes errors, and is readily adapted to automatic plotting methods pro- 
viding computation time is not excessive. 

Figure 8 shows the overall flux-density picture with condensing and diverging pseudo-spherical waves. The depar- 
ture from a Gaussian sphere is dominated by a linear phase gradient of weight (u>m) and a comatic aberration gradient 
of weight (wi 31 ) as described in Appendix B of this document. In the flux-density approach the distorted wavefront per 
se is never identified. 

Figure 9 shows the flux-density picture in the vicinity of the caustic surfaces (co t , oj 2 ) as seen in the (x, z) plane. It is 
particularly important to identify the ray bundle associated with the flux-density mapping and the intercepts of the 
bounding rays of that bundle with the caustic curves. Within the ray bundle there is a region where flux-density values 
superimpose. The introduction of rays outside of the bundle (S/, S r ") enlarges that region. See Ref. 6, p. 259; p. 262. 

When an axial plane wave illuminates a portion of a sphere the results are not as degenerate as for the paraboloid. 
A pseudo-focus appears at (0,0, -c/2), a line caustic (a> 2 ) forms, and a surface caustic (oj 2 ) may be associated with a 
nephroid generating arc rotated about the reflector (z) axis. Figure 10. With paraxial reception, the primary difference 
is a rotation and translation of the caustic surfaces with some distortion (Figure 1 1). It is often instructive to trace the 
progression of the points x(wO» x(a> 2 ) as radial variable (a) is updated for some azimuthal value (f) on the reflector (7). 
A mapping of the flux-density for the sphere with inclined incident plane wave is given as Figure 12. Condensing and 
diverging waves are seen here also, but the fine-structure in the vicinity of the displaced pseudo-focus differs greatly 
from the previous paraboloid example. 

Although this document is primarily concerned with identifying the means for obtaining caustic/flux density infor- 
mation, and verifying some of the formulation, the practical implications for microwaves and optics are also of in- 
terest. The paraboloidal and spherical reflector examples were selected because of the wide application of these sur- 
faces as focussing or collecting devices. Caustic analysis appears to offer an economical means for obtaining high- 
frequency focal region information in three-dimensional space. The analytical generation of coordinates for two 
caustic sheets, with varying degrees of degeneracy, provides insight not readily obtained by the tangent approach 
found in many textbooks and journal articles. Annexation of flux-density greatly enhances the value of the classic 
caustic analysis of wavefronts. Extensions to edge-diffraction, including various approaches such as the geometrical 
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Figure 7. Caustic Surfaces from Flux-Density Singularities (Paraboloid F = 1.0) 
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Figure 9. Caustic Surfaces and Flux-Density (Paraboloid F = 1.0) 
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Figure 10. Caustic Surfaces from Flux-Density Singularities (Sphere C - 1.0) 
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Figure 1 1 . Caustic Surfaces From Flux-Density Singularities (Sphere C = 1.0) 
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theory of diffraction (GTD) and the uniform asymptotic theory of diffraction (UTD), lie beyond the scope of the pre- 
sent discussion. Ref. 24, pp. 67*119. Also, the caustic approach is not placed in competition with other means for ob- 
taining focal-region information (i.e. physical optics for example). An example of such a comparison may be found in 
Ref. 25, p. 477. Also see Ref. 30, p. 266 for the fine structure of a cusp caustic and Ref. 6, p. 263 for examples of 
smooth and cusp caustics. 

The interested reader may wish to review the design of a line-feed built by the Cornell University Center for 
Radiophysics and Space Research for the 1000-foot spherical reflector of the Arecibo Ionospheric Observatory in 
Puerto Rico. Refs. 22, 26, 27, 28. Additional information on a line-source corrector for a large spherical reflector is 
found in Ref. 29. It is noted that every point of the nephroid in Figure 10 is the result of the intersection of two 
rays, and that every point of the line caustic («») is the result of the intersection of a cone of rays. Ref. 23, p. 118. The 
implications of the number of ray intersections at various points in the focal region (cusp of a caustic) with respect to 
field values is discussed in Ref. 31, pp. 1331-1358. Details pertaining to the electromagnetic diffraction in the focal 
region of a spherical mirror under oblique illumination may be found in Refs. 32, 33, and 34. 


CAUSTICS FROM A TENSOR ALGORITHM 


A method of obtaining the equations of caustics by means of a “tensor” algorithm has been found in the literature. 
Ref. 19. The cited source does not contain the derivation of the method. It asserts that 


x, = x, + pX, , (67) 

(«) (7) (7 - «) 

vectorially. Here p is a scalar, X, is a unit vector along a reflected ray, x, is a vector to the reflector 7 , and x, is a vector 
to the caustic (0) of two sheets or branches («,, « 2 ) as illustrated in Figure 13. 

The central problem is to find, explicitly, an equation giving the distance (p) to a caustic (0) that is known in terms 
of a family of tangent rays (S r ) from a reflector ( 7 ). An equation is given in Ref. 19 for the unknown distance 


-gm du a -S2a du a 


p =• 


together with tensors 


J aj3 


and 


g a|3 


: du a 


3Xj 

3 Xi 


3u 0 

3Xi 

3X; 

3u a 

3 Uj3 


G 2 c* du a 


(a = 1 , 2 ) 


(a, P = J, 2) 


(a, (3 = 1, 2) 


( 68 ) 


(69) 


(70) 


In Ref. 19, for the paraboloidal reflector, 

u, = p, u 2 = 6 (71) 

are radial and azimuthal variables of an orthogonal (cylindrical) net, and an inclined plane wave is assumed to be the 
source of rays (S,). 

The reason for introducing the second-rank tensors g a ^ and G a ^ is obscure, and may have been merely a desire 
for a compact expression for distance (p), but this is speculation. A derivation of the present method of obtaining the 
equation of a caustic surface directly from a specification of the source and the reflector, without introducing the no- 
tion of flux-density, is given below. Beginning with equation (67), 


Since X, is a unit vector, 


dx, = dx, + d(pX.) 


X, • dX, = 0 = dx, • dX, 


and 


dXj ) du a = 0 (a = 1,2) 
3u a , 


From this. 


17 \ du a + /p dx j j d V + (Xj)_dp ' 1 _ 

|_\ du a/ \ 9 uq, / duQ,J 

[ 3xj duj + p dx^ duj + ^ + p axjl Tax^u, + axjl 

3uj du 2 3uj du 2 3 u 2 3u 2 J |_ du l du 2 ^ U 2J 


(72) 


(73) 


(74) 


(75) 
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In compact notation, 


(gu + p G,,)/du,\ 2 + (gn + p Gn + gn + pGij)/du\+ (gia + p G 22 ) — 0 (76) 

ydu 2 j ydu 2 y 

If the ratio (dui/du 2 ) is to have only one value, 

(g 2 i + p G 21 + gi 2 + P G 12 ) j - 4 (g„ + p Gu) (g 21 + p G 22 ) = 0. (77) 

For the paraboloid with off-axis plane-wave, 


G 12 = G 21 


(78) 


and 


gu = gu = 0 (79) 

The last equation|is|attributed to the existence of a normal rectilinear congruence in Ref. 19, but appears to be due to 

the selection of an orthogonal (p,0) net. 

Applying both equations, (78) and (79), the simplification obtained for this specialization is 


P 2 (G l2 ) 2 - (gu + pG,i)(g„ + p Gu) = 0 


(80) 


and a quadratic in (p) is recovered: 


P 2 + 


(g 22 Gu + gu G22) P 
(Gu G22 - G.2 2 ) 


+ 


gll §22 


(Gu G22 - G,2 2 ) 


= 0 


The roots of the preceding quadratic, 


(81) 


P = Pt, P 2 , (82) 

are not radii of curvature, but they stretch or compress the unit vector X, (or ,S r ) to locate points of the caustic bran- 
ches («,, « 2 ). In this respect (p) is operationally identical to (R 2 ), associated with the flux-density singularity and 
discussed previously. The derivation, above, circumvents an assumption of Ref. 19 that there exists a du, and, in- 
dependently, a du 2 such that utilization of equation (75) yields the equal lengths (p). Introducing the ratio (du,/du 2 ) 
results in the same quadratic in (p) obtained in Ref. 19, but admits a simultaneous pair dui, du 2 and avoids 
“eliminating” these factors in a separate operation. 

A specific result that utilizes the preceding method is given in Appendix C of this document. The equation for the 
caustic surfaces is obtained for a paraboloid illuminated by an off-axis plane wave. 
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CONCLUSION 


The present document has discussed several methods for determining analytical caustic surfaces. It was shown that 
the mapping of focal regions is an extended study of surfaces including the source or incident wavefront, the reflector, 
and the reflected wavefront in addition to the caustic surfaces, themselves. A review of the classical three-dimensional 
treatment of center surfaces, as applied to distorted wavefronts, was augmented by the flux-density analysis. The lat- 
ter, predicated on the Jacobian of a mapping, obtained the caustic surfaces directly from a specification of the il- 
luminating function and the reflector, to the exclusion of an identifiable physical wavefront. Only the intrinsic Gaus- 
sian, mean, and normal curvature, associated with the incident ray, were required to obtain both caustic branches by 
exploiting a singularity of the flux-density equation. 

This effort progressed until it became apparent that a second document would be needed to display isometric views 
of the caustic surfaces. The computational aspects of the problem have been dealt with using the TURBO PASCAL 
language, but a need for automatic plotting capability with a “hidden-line” feature seemed essential to enhance presen- 
tation, reduce the amount of labor, and guard against human error. An extension of the wavefront approach to a 
Gaussian sphere, rather than a hemisphere, is also being implemented. A discussion of multiple reflection, or higher 
epicycloids, previewed by both the paraboloidal and spherical reflector examples presented here, has been deferred to 
a second document together with such topics as multiple and extended sources. The latter raises the question of 
relative phase,and the phase of the caustic surfaces in general. Speculation regarding “intersecting” and “nearly- 
intersecting” ray density, and other controversial subjects, were also set aside to keep this paper within reasonable 
bounds. 
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APPENDIX A 

Gauss/Seidel Aberrations (Parallel Distortion of Gaussian Wavefront) 


x = o sin f 

z = - (c 2 - a 2 ) Vl 

~ w 020 l 2 

- "ill l 3. 

■ ^040 h 

" w 222 
‘ w 220 ! 1 h 

' w 131 l 2 r 3 

‘ ^311 l 3 *1 


y = - 0 cos f 

(Gaussian Sphere) 

(axial focus) 

(transverse focus) 
(spherical) 

(astigmatism) 

(field curvature/Petzval) 
(coma) 

(distortion) 


where 


Ij = h 2 , I 2 = a 2 , I 3 = h a cos f 
Substituting (A-3) into (A-2), the wavefront z-coordinate is given by 

z = - (c 2 - CT 2 ) 1/2 - (OJ 020 ) o 2 - (C 0 j J J h) 0 cos J - (co 040 ) a 4 - (co 222 h 2 ) a 2 cos : 
- (oj 220 h 2 ) a 2 - (w 131 h) a 3 cos J - (co 3U h 3 ) a cos f 
The wavefront tangents are then obtained as 


X u = [ Sin - cos f, o( c 2 - a 2 y' A 2 - 2(gj 020 ) a - (oo in h) cos f - 4(a> 040 ) a 3 
- 2 (co 2 22 h 2 ) a cos 2 f -2(co 220 h 2 ) a - 3(co 13I h) a 2 cos f - (co 3U h 3 ) cos f | 


and 

x v = [ a cos ?, o sin f, (coj j j h) a sin J + 2(oo 222 h 2 ) a 2 cos J sin f 

+ (w 131 h) a 3 sin f + (co 311 h 3 ) a sin J ] 

The second partial derivatives are 

\m = ( ° 2 ( c2 " ° 2 y 312 + (c 2 * ^y' h - 2 H20) ~ 12(co 040 ) a 2 

-2(co 222 h 2 ) cos 2 f - 2 (co 220 h 2 ) -6(oo ]31 h) a cos ? ] 
x uv = I cos sin ( w m h) sin f + 4(^222 0 sin f cos f + 3 (cOjji h) a 2 sin f 

+ (co 311 h 3 ) sin f ] 

x vv = [ - a sin f, a cos f, (a> m h) a cos f + 2 ( CJ 2 22 h b ° 2 ( cos2 ? “ sin 2 f) 

+ (w 131 h) a 3 cos f + (co 311 h 3 ) a cos f ] 
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APPENDIX B 

Gauss/Seidel Aberrations (Normal Distortion of Gaussian Wavefront) 


Ns = - [ a sin J , - o cos f, - (c 2 - a 2 ) 1 ' 4 ] / c 

is taken as the normal to the Gaussian sphere of radius (c). The wavefront coordinates are given by 

x = o sin f + (^020^ ° 3 sin f + ( w Hl h ) ° 2 sin f cos f + (o; 040 ) ^ sin ? 

c c c 

+ (CO222 h 2 ) V* sin f cos 2 f + (CO220 h 2 ) a 3 sin J + (a) 131 h) a 4 sin f cos f 
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The wavefront tangents are then obtained as 
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APPENDIX C 

The Caustic Surfaces of a Paraboloid and Inclined Plane Wave 


A complete verification of the caustic equation obtained in Ref. 19 for the paraboloid and inclined plane wave is 
tedious. Some of the important factors are given here for convenience, together with the final result. In Ref. 19, the 
parametric equations of the reflector surface are given as 


Xj = — (p 2 - 4a 2 ) 

1 4a 


(C-l) 


x 2 = p cos 6 


(C-2) 


x-j = p sin 6 (C-3) 

The details pertaining to the reflected ray are omitted entirely here since they are easily obtained from 

S r = S. - 2(N • Sj) N (C-4) 

The components of a reflected unit vector are 
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4a 2 + p 2 


J 


- 4a p sin 6 cos <p + 2 p 2 sin 2 6 sin <p \ /\ 


4a 2 + p 2 


From the preceding Xj, the partial derivatives 

H(X CUS ^ \J 

9uj 9 p 


/\ 


9X i 9X i 4a cos <p (p 2 sin 8 tan <p - 4a p - 4a 2 sin 0 tan <p) i + 


(4a 2 + p 2 ) 2 


(C-5) 


/V 


4a cos (p 2 cos 8 + 4a p sin 0 cos 8 tan - 4a z cos 0) j + 

(4a 2 + p 2 ) 2 


/\ 


4a cos <p (p~ sin 8 + 4a p sin z 8 tan \p - 4a z sin 0) k 
(4a 2 + p 2 ) 2 


(C-6) 


and 
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3X i 3X i -4a p cos 0 sin <p /\ 

^ 4a 2 + p 2 

p cos (4a sin 9 + 2 p cos 2 0 tan <p) ^ + 

4a 2 + p 2 

p cos <p (- 4a cos 9 + 2 p sin 2 0 tan <p) it" + 
4a 2 + p 2 


give 

Gj j = 16a 2 (cos 2 <p + sin 2 0 sin 2 <p) / (4a 2 + p 2 ) 2 

Gj 2 = G 2 ] = 8a p cos 9 sin <p (p cos <p + 2a sin 0 sin <p) / (4a 2 + p 2 ) 2 

G 22 “ 4p 2 [4a 2 + p 2 - (2a sin 0 sin <p - p cos <p) 2 ] / (4a 2 + p 2 ) 2 

The reflector tangents 


3 ^ 

3u t 


p y\ 

— L = — i + cos 0 j + sin 0 k 

3p 2a 


3x- 3x- /\ 

— - = — - = - p sin 0 j + p cos 0 k 

3u 2 30 

together with the partial derivations of X i? give 

gll = 2(p sin 0 sin <p - 2a cos <p) / (4a 2 + p 2 ) 


812 g 21 


= 0 


g 22 = 2p 2 (p sin $ sin <p - 2a cos <p) / (4a 2 + p 2 ) 

The equation for the coordinates of the caustic surfaces of the paraboloid with incident plane wave inclined 
angle (<p) is obtained as 

(2a + X|) cos 


= x. 


+ 


1 - 


(2a cos <p - si 


)s <p 

— — - Q 

sin <p) 


i 


r x 2 " 


x 2 - — 

Q j + 

I 2 2a 


x 3 “ 


x 3 (2a + Xj) sin <p 

2a (2a cos <p - x 3 sin <p) 


I 


/\ 

k 


where 


0 = 2a + Xj - cos <p( x ] cos <p + x^ sin <p) ± (2a + Xj) sin <p 

and the point (Xj, x 2 , x^) lies on the paraboloid. 


V I _ / X] cos <p + x 3 sin ip V 

\ 2a + Xl / 


(C-7) 

(C-8) 

(C-9) 

(C-10) 

(c-n) 

(C-12) 

(C-13) 
(C- 1 4) 
(C- 15) 
at an 


(C-16) 

1 

(C-17) 
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